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Abstract. We present an algorithm for the construction of conformal coordinates in the interior of 
a spherically symmetric, collapsing matter cloud in general relativity. This algorithm is based on the 
numerical integration of radial null geodesies. As an application we generate conformal diagrams 
for collapsing spherical dust clouds and analyze the causal structure of the resulting spacetimes. 
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INTRODUCTION 

General relativity predicts that sufficiently massive, cold fluid stars cannot exist in hydro- 
static equilibrium and, hence, such stars must undergo complete gravitational collapse 
(see [1,2] and references therein). If a trapped surface forms during the collapse, then 
the singularity theorems [1] also imply that the final state includes a spacetime singular- 
ity. One of the interesting questions is whether this singularity is censured, that is, if it is 
hidden inside a black hole, or if it is naked in the sense that light rays emanating from the 
singularity can reach observers which are arbitrary far from the collapsing star. Accord- 
ing to the weak cosmic censorship conjecture [3] the former case must occur. A stronger 
form of this conjecture [4] states that the future Cauchy development of generic, regular, 
asymptotically flat initial data gives rise to a globally hyperbolic spacetime, implying 
that the singularities in this spacetime are not visible to any observer. 

On the other hand, it is known that the collapse of a spherical dust cloud leads to the 
formation of shell-focusing singularities which may be visible to local observers, and 
even to distant observers [5, 6, 7]. However, it is not clear whether or not these examples 
provide a serious threat to the cosmic censorship conjecture. Indeed, it remains to be 
seen if the nature of the singularity is unchanged under small perturbations of the initial 
data and the equation of state. For reviews and a recent book on this subject we refer the 
reader to [8, 9] and [10], respectively. 

In this article, we describe a numerical method for the construction of conformal 
diagrams in the spherical symmetric dust collapse. We apply this method to different 
initial data and analyze whether the resulting singularity is censored or naked. Our 
method should be generalizable to the spherical collapse of more general matter and 
is expected to be useful as a tool for an understanding of the dependency of the nature 
of the singularity (censored or naked) on the initial data. 
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THE MODEL 

As a simple collapse model we consider a self-gravitating, spherically symmetric dust 
cloud. In terms of co-moving, synchronous coordinates [11] this model is de- 

scribed by the one-dimensional mechanical system 

^f{T,Rf + V{r{T,R),R)^E{R), V{r,R) := (1) 

where R labels the different dust shells and where the function t i-> r{x,R) describes the 
evolution of the areal radius along dust shell 7? as a function of proper time. It is assumed 
that r{0,R) = R initally; that is, R coincides with the areal radius of the shell R at the 
initial time T = 0. The energy function E{R) = vo{R)^/2 — m{R) /R is determined by the 
initial velocity and density profiles vq{R) := r{0,R) and po(i?), respectively, where m{R) 
is the Misner-Sharp mass function [12] obtained from 

R 

m{R) = 4nG J po{s)s^ds, (2) 


with Newton's constant G. Once the function r{x,R) has been determined, the time- 
evolution of the metric g, the four- velocity u and the density p are obtained according 
to the following formulae: 

g = -dT^ + ^^^^dR^ + r{T,Rf{d^^ + .\nHd<p\ 

d . ^ , R \^ I 

u = ^, pM=poW(^j 

Here and in the following, a dot and a prime denote partial differentiation with respect 
to T and R, respectively. 

Eq. (1) can be solved exactly [7]. Since we are interested in the collapse of a dust 
cloud with positive density whose initial data does not contain any apparent horizons, 
we make the following assumptions 

P>0, vo<0, -^<E<0. (5) 
In this case, the solution to (1) is given by 

r{T,R) = ^F-i (^F{p{R)) + ^c{R)p{R)H^ , 

where p{R) := E{R)/V{R,R) E (0, 1] is the ratio between the total and initial potential 
energies, c{R) :— 2m{R)/R^ is proportional to the mean density within the dust shell 
R and where the strictly decreasing function F is given by F : [0,1] — )■ [0,;r/2],x i-> 
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(3) 
(4) 



y(^(T^^^ + arccos In the following, we also assume that the functions p and c are 
subject to the inequalities 

p' > 0, c' < 0. (6) 

These conditions are fulfilled, for instance, if the density profile is decreasing and the 
initial velocity is zero. They imply the absence of shell-crossing singularities [7]. There 
remains a shell-focusing^ singularity at points where r(T,i?)/i? = 0, in which case 



T=T,(i?) = 



f-nP(^)) 



This corresponds to a curvature singularity since the Ricci scalar diverges as a conse- 
quence of Einstein's field equations and the divergence of p at those points. 

In the next section we describe a method for generating conformal diagrams for the 
spacetime metric (3) subject to the assumptions (5,6). By definition, in- and outgoing 
radial null geodesies in such a diagram are given by straight lines with slope — 1 and 
1, respectively, which nicely displays the causal structure of the underlying spacetime. 
In particular, we are interested in determining whether the shell-focusing singularity is 
censored or naked. 



CONFORMAL DIAGRAMS 

In this section we describe a method for constructing conformal coordinates {T,X) 
for the spacetime described by Eq. (3). We assume that the cloud has a finite initial 
radius Ri > 0. For R> Ri the density is zero and according to Birkhoff's theorem the 
spacetime metric (3) must reduce to the Schwarzschild solution whose causal structure 
is known. Therefore, it is sufficient to construct conformal coordinates for the region 
Q. := {(t,/?) :0 <R< Ri,0 < T < inside the dust cloud. Such coordinates can 

be obtained by first solving the advection equations 

u = -ru', v = yv', r-=^ — , U) 

r 

with appropriate initial and boundary data for the null coordinates U and V and then 

setting T = {V + U) /2, X = {V — U)/2. Eqs. (7) in turn can be reduced to ordinary dif- 
ferential equations by the method of characteristics: the functions V and U are constant 
along the integral curves belonging to the null vector fields 

-~ dz ^dR' -'dz^^dR' 

respectively. These integral curves describe in- and outgoing radial null geodesies. Since 
Y > on Q., the integral curves belonging to k^'"^ cannot cross at points away from 
the singularity. The same statement holds for the integral curves to k^""'\ Furthermore, 



See [13] for the distinct physical properties of shell-crossing and sheU-focusing singularities. 
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a careful argument [6, 7] shows that under appropriate regularity assumptions on the 
initial data po and vq, the ingoing radial null geodesies cannot cross on the singularity 
y {{t,R) : < R < Ri : T = Ts{R)}, and that the outgoing radial null geodesies 
cannot cross on y except at the point — (0, Tj(0)). Therefore, we can construct U 
and V by specifying their values on the singularity y, by specifying V at the center 
line ^:={(T,i?) : R = 0,0 < T < z^} and by specifying U at the surface of the cloud, 
^ := {{t,R) : R = Ri,0 < T < T,(7?'i)}. 

We first apply our method to the case of a homogeneous dust cloud with zero initial 
velocity and then describe a numerical algorithm which is applied to the generic case. 



Special case of a homogeneous cloud with zero initial velocity 

For a dust cloud with homogeneous density and zero initial velocity, c{R) = cq 
is constant and p = I on Q.. The assumptions (5) require < cq < 1. In this case, 
Y— a/ 1 — cqR'^ / {^/cqt) and the coordinate transformation 

^/coT = F ^cos^ , y/coR = sm{X), < 7 < ;r,0 < Z < arcsin(^i?i) 

gives 

cos2 (I) \dT dxj ' " cos2 (I) {dT^dXj- 

Therefore, a solution to the advection equations (7) is U = T + X and V = T —X. The 
corresponding conformal diagram for the choice cq — 0.75 is displayed in Figure 1. 
Here, the apparent horizon is the boundary between the points for which the areal radius 
r increases along outgoing radial null geodesies, k^'^'*'^ [r] > 0, and the points for which 
j^{out) ^ 0. A small calculation reveals that this surface is determined by the simple 
equation r(T,i?) = 2m{R). For the homogenous case with zero initial velocity considered 
here this reduces toT — K — 2X. Since outside the cloud the apparent horizon coincides 
with the event horizon, the event horizon inside the cloud is determined by the outgoing 
null ray passing through the intersection of the apparent horizon with the surface of the 
cloud. 



The generic case 

For the generic case, we apply the algorithm described above in order to determine 
the functions U and V on the inside region Q.. To this purpose, we select a point p 
inside Q.. Then, we numerically integrate the future-directed in- and outgoing radial null 
rays emanating from p. The ingoing ray ends at either the center of the cloud ^ or at the 
singularity y. The outgoing ray ends at y or at the surface of the cloud, ^. From these 
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FIGURE 1. Conformal diagram for a homogenous dust cloud with zero initial velocity, cq = 0.75 and 
Ri = l. The Unes denoted by "AH" and "EH" refer to the apparent and event horizons, respectively. 



two end points we define the coordinates {U,V) associated to the point p according to 



u 



5^ 



1- 



R 
Ra 



Notice that this choice implies that the singularity ^ \ T^^ is determined by the line 7=1 
and < Z = RfRi < 1 in the conformal diagram. On the other hand, it turns out there 
are infinitely many outgoing null rays emanating from T^^. These null rays have the same 
value of y = 1 and different values for U . Therefore, the point unfolds to a portion 
of the straight line T = \ —X m the conformal diagram. Two examples are shown in 
Figures 2 and 3. In the first case, the singularity is hidden inside the black hole while in 
the second case a portion of is visible from future null infinity. The results are based 
on the following choice for the functions c and p\ 



c{R) = CO 



^ 6 fRY 3,^ ^JR 



p{R)^Po + {pi-po)^^, (8) 



where the parameters cq, a, po and pi are subject to the inequalities 

7 



0<a<l, 0<co< 



6fl + 4' 



0<po<Pi<l, 
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FIGURE 2. Conformal diagram for the model described in Eq. (8) with the parameter choice cq = 0.65, 
a = 0.7, po=0.\ and pi = 0.9. The lines denoted by "AH", "EH" and "CH" refer to the apparent, event 
and Cauchy horizons, respectively. In this case, the singularity is hidden inside the black hole region 
since all light rays emanating from it end at the singularity T = 1. 



which guarantee the satisfaction of the assumptions (5,6). Furthermore, c{R) is such that 
the density p converges to zero as R ^ R\ approaches the surface of the cloud. 

For the numerical integration of the light rays we first rescale the vector fields ^('"^ 
and k^'"*'^ by the factor (1 + y2)~V2 and then solve the system of ordinary differential 
equation 

^ 1 dR ^ r{^,R) 

dX ^\ + y{x,RY' dX ^l + 7(T,i?)2' 

This rescaling corresponds to a reparametrization of the light curves and offers the 
advantage of avoiding stiff coefficients. Therefore, adaptive or implicit time integration 
is not necessary. We numerically integrate Eq. (9) using a fourth-order Runge-Kutta time 
integrator (see, for instance, [14]) with fixed step size h. The results for the apparent 
horizon corresponding to different values for h are shown in Figure 4, and indicate 
convergence of our method. We have also verified the convergence of the earliest point 
of the image of in the conformal diagram, from which the Cauchy horizon emanates. 
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FIGURE 3. Conformal diagram for the model described in Eq. (8) with the same parameter choice as 
in the previous figure except that cq =0.45. In this case, there exists hght rays emanating from Tj which 
arrive at the surface of the cloud earlier than the apparent horizon. Therefore, a portion of this singularity 
is visible to distant observers outside the black hole region. 



CONCLUSIONS 

We presented a simple numerical method for determining the conformal diagram inside 
a collapsing, spherical dust cloud with generic initial density and velocity profiles. This 

provides a valuable tool for studying under which conditions the singularity is hidden 
inside a black hole and under which conditions it is visible from future null ifninity. Our 
method should also be applicable to spacetimes describing the spherically symmetric 
collapse of clouds with nonzero pressure since it is based on the integration of radial null 
rays. However, unlike to present dust case, the spacetime metric for generic initial data 
can hardly be obtained in analytic form for realistic equations of state, and a numerical 
integration is probably necessary for obtaining the solution to the field equations. 
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FIGURE 4. The location of the apparent horizon in the conformal diagram for different resolutions. 
The corresponding convergence factor oscillates about a value close to three, indicating al least first order 
convergence. The initial density and velocity profiles are the same as in the previous figure. 
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